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With the upgrade of current gravitational wave detectors, the first detection of gravitational wave 
signals is expected to occur in the next decade. Low-latency gravitational wave triggers will be 
necessary to make fast follow-up electromagnetic observations of events related to their source, e.g., 
prompt optical emission associated with short gamma-ray bursts. In this paper we present a new 
time-domain low-latency algorithm for identifying the presence of gravitational waves produced by 
compact binary coalescence events in noisy detector data. Our method calculates the signal to noise 
ratio from the summation of a bank of parallel infinite impulse response (IIR) filters. We show that 
our summed parallel infinite impulse response (SPUR) method can retrieve the signal to noise ratio 
to greater than 99% of that produced from the optimal matched filter. We emphasise the benefits 
of the SPUR method for advanced detectors, which will require larger template banks. 
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I. INTRODUCTION 

The interferometric gravitational wave (GW) detectors 
LIGO pQ, and Virgo [2] have reached a sensitivity at 
which the detection of GWs is possible. The LIGO de- 
tectors are currently undergoing a major upgrade to Ad- 
vanced LIGO, for which the sensitivity will be improved 
ten fold relative to Initial LIGO [3]. Hence Advanced 
LIGO will be able to detect GW (GW) sources within a 
volume of space one thousand times larger than that of 
initial LIGO, out to -200-300 Mpc @]. 

The emission of GWs produced by compact binary co- 
alescence (CBC) can be modelled with a high degree ac- 
curacy [5]. When two compact bodies, such as neutron 
stars or black holes are in orbit, Einstein's equations pre- 
dict the generation of GWs. As the bodies spiral towards 
each other a GW is created that increases in frequency 
over time until the bodies merge, following what is known 
as the inspiral waveform. Ground based detectors have 
frequency passbands that allow them to be sensitive to 
the final stages of such events up to a total system masses 
of several hundred M . 

Neutron star binary mergers are widely thought to be 
the progenitors of short hard gamma-ray bursts (short 
GRBs) [6j [7] . The delay between the final GW emission 
and the onset of the GRB is estimated to be as short 
as 0.1 seconds or as long as tens to hundreds of seconds 
[H [H] . The electromagnetic emission of the GRB event is 
not well understood. Related to the initial GRB there is 
thought to be a prompt emission in the X-ray and optical 
wavelengths followed by a delayed afterglow of cascading 
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wavelengths. Prompt optical emission may occur tens 
to hundreds of seconds after the initial burst. The low- 
latency detection of the GW associated with a neutron 
star merger could lead to the localisation of a GRB source 
event on the sky, enabling fast moving telescopes to ob- 
serve the prompt optical emission. Data collected from 
a multitude of sources — GWs, gamma-rays, X-rays and 
optical counterparts of the GRB — will lead to maximum 
insight into these highly energetic events. 

The standard strategy for searching for the existence 
of inspiral waveforms in the detector data is based on 
matched filtering [5] (and references therein). This 
method, based on Wiener optimal filtering, is a correla- 
tion of an expected inspiral waveform template and the 
detector data, weighted by the inverse noise spectral den- 
sity of the detector [TU]. In order to save computational 
costs, this correlation is performed in the frequency do- 
main, via a Fourier transform of a finite segment of detec- 
tor data. In previous LIGO searches, the detector data 
is split up into "science blocks" , which are further di- 
vided into "data segments" chosen to be at least twice 
the length of the longest waveform in the template bank 
|11) . Each proceeding data segment is chosen to over- 
lap the previous one by 50%. Each segment therefore 
must be matched filtered in a time that is half the length 
of the segment for a real-time analysis, that is, the fil- 
ter output rate is equal to the data input rate. In this 
case, the matched filter process has a minimum latency 
(from signal arrival to signal detection) that is propor- 
tional to the longest template (see [T^] for more details) . 
Advanced LIGO will have an increased bandwidth over 
Initial LIGO, with the lower bound dropping from 40 Hz 
to 10 Hz [3]. GW signals from CBC events spend much 
more time at these lower frequencies. Hence waveforms 
used for matched filtering in Advanced LIGO will be 
much longer (1000s of seconds). This in turn means the 
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segment length will be increased, further increasing the 
latency. The latency of this method to produce GW trig- 
gers is longer than the time to onset of prompt optical 
emission after coalescence (10s to 100s of seconds). After 
this amount of time, the early electromagnetic counter- 
part of a GRB event will be significantly faded, and may 
be missed by telescopes altogether. 

A low-latency GW detection method is required to 
trigger follow-up electromagnetic observations of the 
prompt optical emission. So far two frequency domain 
methods have been developed to solve this issue. The 
VIRGO group has produced a low-latency pipeline based 
on Multi-Band Template Analysis (MBTA) |13j . and 
LIGO is also working on a new method, Low-Latency On- 
line Inspiral Data analysis (LLOID) method. In MBTA 
the matched filtering technique is split over two frequency 
bands, and the output is coherently added, reducing la- 
tency. A latency of less than 3 minutes until the availabil- 
ity of a trigger using this method has been achieved [12] • 
Low-latency in the LLOID method is achieved by first 
down-sampling the incoming data into multiple streams 
and then applying frequency domain finite impulse re- 
sponse (FIR) filters [14 . The computational cost of this 
pipeline is reduced by decreasing the number of templates 
via singular value decomposition |15) . 

We introduce a new method to detect CBC signals in 
the time domain using infinite impulse response (IIR) 
filters. Approximating an inspiral waveform by a sum- 
mation of time shifted exponentially increasing sinusoids 
enables us to construct a bank of parallel single-pole IIR 
filters. Each IIR filter acts as a narrow bandpass filter. 
When each appropriately delayed IIR filter is added the 
coherent output approximates the matched filter output 
of the exact waveforms. We call this the summed par- 
allel infinite impulse response (SPUR) method. Figure 
[T] visually demonstrates the idea of using a bank of IIR 
filters as narrow bandpass filters. For a full explanation 




FIG. 1: A schematic overview of the SPUR method. The in- 
put is split into different channels, time delayed by an amount 
d, then passed through a narrow bandpass IIR filters, each 
with a different central frequency /. Finally the output of 
each individual IIR filter is summed, giving the output of the 
SPUR method. 

of the mathematical principles, see [12] . In this follow up 
paper, we numerically address the issues essential to the 
practical use of this method for the upcoming advanced 



detectors. We calculate the filter coefficients and demon- 
strate via numerical simulations how well our method ap- 
proximates the optimal matched filter as a function the 
number of filters per bank using a range of parameters. 
We also show that the detection rate of the SPUR method 
is very similar to that of the matched filter method. It 
has been shown theoretically that in order to get the 
same latency as the SPUR method, the frequency do- 
main matched filter method would require greater com- 
putational resources [12) . 

The structure of this paper is as follows: In section 
|TT| we will go through the formal introduction of the in- 
spiral waveform and matched filtering, and how to get 
from the continuous frequency domain matched filter to 
the time domain discrete matched filter. This will lead 
to a demonstration on how it is possible to approximate 
an inspiral signal by a sum of exponentially increasing 
sinusoids. The methodology is explained in Section |III| 
and will cover how we set up our simulation to test the 
efficiency of the SPUR method as opposed to the fre- 
quency domain matched filter. Section IV will analyse 



the results of the simulation and SectionTvl will discuss 
the implications of these results for advanced detectors. 



II. METHODOLOGY 

Gravitational wave interferometers output the strain 
induced by gravitational waves incident on the detector, 
as well as inherent noise. In unitless strain, the detector 
output will be, 



n(t) 



if signal is absent 



I n[t) + h(t) if signal is present 

where n(t) is the noise inherent in the detector. The 
sensitivity of the instrument can be characterized by the 
(one-sided) strain power spectral density S n (f), 
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(n(f)n*(f')) = -S n (f)6(f-f) 



(2) 



where the tilde represents the forward Fourier transform, 

poo 

q(f) = / q(t)e- 2 ^dt. (3) 



A. The Inspiral Waveform 

The gravitational-wave strain incident at the interfer- 
ometer is given by 

h{t) = F+(0, 0, $)h+ (t) + F x (9, 0, i/,)h x (t) (4) 

where the detector response functions F+ and F x are 
functions of (6, cj>) - the standard spherical polar coordi- 
nates measured with respect to the Earth's fixed frame, 
and if) is the polarisation angle. The detector response 
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function can be found in [T6] . The + and x polarisations 
of the waveform are, 



h+(t) = 



1 



A(t) cos <j>(t) 



h x (t) = (cos l) A(t) sin (j>{t) 



(5) 
(6) 



For non-spinning binaries with a chirp mass 
M = ((mim 2 ) 3 /(mi + rf"^)) 1 / 5 in the range of 1 — 3M 
- we will hereafter assume — the waveforms can be 
modelled to very high accuracy using the Restricted 
post-Newtonian (PN) expansion fTTHE)] in the LIGO 
band (assumed to be 10-1500 Hz for advanced LIGO). 
For restricted waveforms, only the leading order of the 
amplitude A(t) is taken, 



A{t) 



GM 



(-. 



t r - t 



-1/4 



Dc? \5GM/c 3 / 
and the post-Newtonian phase <p(t) is given by 



(7) 



t c - t 
GM/c 3 



5/8 



+ higher order terms (8) 



In addition to the source masses mi, ma, there are sev- 
eral unknown parameters; the time of coalescence t c , the 
phase at coalescence <po, distance from observer to source 
D, the inclination angle of the binary's orbital plane rel- 
ative the line of sight l, and the polarisation angle ip. 
However by using the linear combination trigonometric 
identity, one can re-express the strain Q by splitting 
the scaling factor due to distance, sky location and ori- 
entation to the mass dependant time evolution of the 
waveform 1201. 



h(t) = 1 " f f PC [h c (t) cos 4> + h s (t)sin(j) } 

where the scalar factor D e g is, 

D 



D, 



eff 



Fl (1 + cos 2 t) /4 + F 2 (cost) 



which gives (f>o, an unknown phase as, 



o = 



arctan 



F x (2 cost) 



COS^ I 



(9) 



(10) 



(11) 



We now define the terms h c and h s as the waveform at 
4>o = and (/> w /2, scaled at 1 Mpc as the so called "cosine" 



and "sine" phases 



h c (t) = A 1Mpc (t) cos (j){t) 
h s (t) = AiM P c(*)sin0(t) 

B. The Matched Filter 



(12) 
(13) 



The matched filter Q is a linear operator that max- 
imises the ratio of "signal" to "noise" present in the de- 



tector data s 



z{t) = 2 



It is denoted by, 



~<f)Q*U) 

Sn(\f\) ' 



,2-niftjf _ 



df = (s(t) \Q) (14) 



Where we have also defined the inner product (a \b) . The 
signal to noise ratio (SNR) is generally defined as the 
ratio of observed filter output to it's expected root-mean 
square flucations or standard deviation, 



SNR = 



(15) 



Note that in the absence of a signal, (SNR) = and 
((SNR) 2 ) = 1 independent of the normalisation of the 
filter Q. 



C. Two-Phase Filter 

A convenient way to search for the unknown phase 
constant O is to filter both phases h c and h s separately 
and then combined to form a complex signal. The two- 
phase filter is defined as, 



z{t) = (s(t)\h c ) +i(s{t)\h s ) 

_ , Hf)K(f) 
s n (f) 



df 



i2 r s(f)ht(f) Mft 

J-oo S n (\f\) 1 



(16a) 
(16b) 



The advantage of using the phases h CtS is that in the sta- 
tionary phase approximation [23], h c and h s are exactly 
orthogonal ((h c \h s ) = (h s \h s ), (h c \h s ) = 0). It then 
follows, h c (f) — ih s (f) for / > 0. Generally, this is 
applied to (161 to give the two-phase matched filter as, 



z(t) = 4 



Sn(\f\) 5 



(17) 



However in this paper, we prefer to maintain the form of 



the two-phase filter in (16 1. In convention with the field, 



the amplitude signal to noise ratio of the (quadrature) 
matched filter is defined as the absolute value of the two- 
phase filter, divided by a normalisation constant that is 
equal to standard deviation of the real and imaginary 
parts of the two-phase filter, 



Pit) 



where a 2 is, 



a z = 2 



k(f) 



df = (h c \h c ) 



(18) 



(19) 



Note that in the in the absence of a signal (just noise), 
the SNR p ( 18 1 is Rayleigh distributed with mean 
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and variance 1, which is identical to the Chi-distribution 
with two degrees of freedom (one for each of the phases). 
This of course implies that the SNR squared, p 2 is Chi- 
square distributed with two degrees of freedom. Hence 
the probability of finding an SNR value greater than p* 

is eh, 



p{p 2 >pD 



(20) 



D. Digital Time Domain Filtering 



The two-phase matched filter 16 is a cross correlation 
of phase h cs (t) and the detector output s(t), weighted 
by the inverse noise spectral density S n (f). By defining 
the quantity x as the over-whitened strain data, 



x(t) = 



~<f) 
S(f)' 



,2mft 



df 



(21) 



we can use the cross-correlation theorem to define the 
two-phase matched filter in the time domain, 

z(t) = 2 [ x(t')h c (t' - t)d£ + i2 [ x(t')h s (t' - t)d£ 



= 2 



x(t')h(t' - t)dt' 



(22) 
(23) 



where h = h c {t) + ih s (t) = A{t)e^ . 

The discrete form of the continuous time domain 
matched filter (1231) is, 



(24) 



z k = 2 ^2 Xjhj-k&t 

j=—oo 



where t — kAt. In practise, the inspiral waveform tem- 
plate hi is bounded (because the detector is only sensitive 
over a bandwidth), and the summation becomes finite, 
making this a finite impulse response (FIR) filter. 



E. Infinite Impulse Response Filter 

Now let us introduce an alternative digital filter, the in- 
finite impulse response (IIR) filter. The difference equa- 
tion of a general IIR filter is, 



Vk 



N 

E 



a„yk-n 



M 



(25) 



outputs are fed back into the filter. However digital IIR 
filter design is a more complex process than FIR design. 
Obtaining the coefficients is usually done by first con- 
structing an equivalent analog filter and applying well- 
known methods, such as the bi-linear transform or im- 
pulse invariance. Multiple IIR filters used together have 
different forms, such as direct form I & II, cascade (se- 
ries) and parallel. In a series configuration, the overall 
transfer function is the multiplication of each IIR filter 
transfer function. In a parallel bank of IIR filters, where 
the output is summed together, the overall transfer func- 
tion is the summation of the different transfer functions. 

First, let's analyse the simplest single-pole IIR filter. 
The difference equation of this filter is 

Vk = aiVk-i + b Q x k . (26) 
A solution to this first-order linear inhomogeneous dif- 




FIG. 2: A signal processing schematic showing the flow of 
data through a digital single-pole IIR filter. The input, x k 
is multiplied by a complex constant i>o, then added to the 
previous output that has been multiplied by another complex 
constant ai, resulting in the current output y k . It should be 
noted that this filter, in principle, should be have been run 
forever. 



ference equation is 



Vk 



k 

E 

3=— oo 



; k—j 

Xjboa 1 . 



By defining the complex coefficient a\ in the form, 



(27) 



(28) 



and comparing (24) and (|27[), it is easy to see that the 



output of the simple filter|26| is the cross-correlation 
of x k and complex sinusoid u n with frequency us and a 
magnitude that increases with an exponential factor 7 
for n < 0: 



u n = 6 e (7+Iu) ™ At e(-n) 
where 0(— n) is the Heaviside function. 



(29) 



where y k is the filter output at time step k, (t = kAt), x k 
is the filter input, and a's and 6's are complex coefficients. 

Examples of IIR filters in common usage are Cheby- 
shev, Butterworth and elliptic filters. IIR filters use much 
less computational resources than an equivalent FIR fil- 
ter. This is because they have "memory" — the previous 



F. Approximation to an inspiral waveform 



Since cj)(t) is not linear in time, a complex sinusoid ( 29 ) 
cannot approximate the h c s phases of the inspiral wave- 
form h{t) — A{t)e 1 ^^ . However we can easily linearise 



the phases by a first-order Taylor expansion about the 
time tf: 



A{t)e 



i<t>(t) 



A(t?)e 



i0(tn+^(tD(t-t(*). 



(30) 



since the amplitude A(t) does not increase at the same 
rate as </>(£), only a linear expansion of <p{t) is required. 
Multiplying by the window function e 7,( ^ t_t, - ) 9(t; — t) 
makes this approximation an exponentially increasing 
constant frequency complex sinusoid with cutoff time tf. 

Ul {t) = A(^)e i ^( i *) + *cn(*'-*r)) e (7 ( +^(tr)(*-* i )e(^ _ t y 

(31) 

The expansion point t ; * is chosen to be near the cutoff 
time, i ; * = ti — aTi, where a is a tunable parameter and 
the interval T} is the duration in which the approximation 
is valid: 



1 



\-4>{ti)T(\ =e<l 



(32) 



and e is a tunable parameter chosen to be to small. Equa- 



tion (|31|) implies that the coefficient bo for the Zth com- 

A ( t *y(4>(t;)+4>(t;)(ti-tD) (33) 



plex sinusoid is, 
bo,i 



and the frequency uji = 4>{t*). 

In this paper, we chose the cutoff time t\ of the first 
sinusoid to correspond to the time at which the waveform 
has the highest frequency detectable by the LIGO detec- 
tor band. The next sinusoid is chosen by moving to an 
earlier time, = t\ — T\. Since we want the Zth sinu- 
soid to be mostly present on the interval ti — Ti < t < ti, 
we choose the damping factor to be 7/ = f3/Ti, where /3 
is a tunable parameter. This procedure is repeated until 
the time t\ corresponds to a time in the waveform that 
has frequency below the LIGO detector band. Hence the 
number of sinusoids is dependent on the value of e, the 
rate of frequency change ((/>(£)), which is dependent on 
the masses of the system, and the detector bandwidth. 
For more information on this procedure, see |12j . 

We can now approximate the phases hit) = A{t)e l ^ < ^ > 
by an addition of a series of damped sinusoids u(t) with 
cutoff times ti : 

A{t)e^ (t) ~U(t) = J2 u i(t) 
1 

= ^ 6o,ie (7,+iWi)(i - tt) e(t i - t). (34) 

Figure [3] shows an illustration of how damped constant 
frequency sinusoids can add to give an inspiral like wave- 
form. 




FIG. 3: An illustrative diagram demonstrating the ability 
to linearly sum exponentially increasing constant frequency 
sinusoids to approximate an inspiral like waveform. The top 
three panels (a-c) show three example sinusoids with different 
damping, frequency and cutoff time factors. Panel (d) shows 
the linear addition of all the sinusoids (at different scales). 
Panel (e) shows the exact inspiral-like waveform. Note that 
this figure is only for illustrative purposes. 
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Here the cutoff time is incorporated by running 
1 filter on a delay, di — t[/At. The output of the Ith 



eac 

filter at time k is 



Uk,l = aiAUk-1,1 + b ,iXk-di 



(35) 



The linear summation of the output of all filters is the 
cross-correlation of the data x and the approximate wave- 
form U(t) in p4l): 



Zk 



2AtJ2 



(36) 



Here z is equivalent to the value computed by the discrete 
time domain two phase filter ( |24[ ) when using a template 
h(t) = U(t). From equation ( |18[), i t follows that the ab- 
solute value of the summation p6| divided by ajj is the 
SNR, which we term the output of the Summed Paral- 
lel Infinite Impulse Response (SPUR). The normalisation 
factor (Ju is defined as 



= 4 



SnU) 



df. 



(37) 



G. Summed Parallel IIR filtering 



Each complex sinusoid ui(t) in equation (34) can be 



searched for in the data x using the single pole IIR filter 



Where Usa(f) is the Fourier transform of the real part of 
U(t), (which approximates h c {t)). The similarity of the 
SPUR output and the matched filter output will depend 
on how well U(t) approximates the given template. 



G 



III. IMPLEMENTATION FOR PERFORMANCE 
TESTING 

A. IIR bank construction 



somewhere within those 4 seconds. At a sample rate 
of 4096 Hz, the analysis period is iV ana i ys i s = 16834 data 
points long. In our simulation, we find d max — 4081683 
and 27„^ x = 149432, resulting in A input = 4247499. 



To confirm the ability of the SPUR method to recover 
a good SNR, it is first required to show that the approx- 
imate inspiral waveform (34) is a good "match" to the 
theoretical inspiral waveform Q. We define the overlap 
A as the inner product of the approximate waveform U 
and the template h: 



A = 



u 



(h\U) 



(38) 



We initially approximate a canonical 2PN 1.4-1.4 Mq 
inspiral waveform band limited to 10-1500 Hz using the 
value of the tunable parameters e, a and f3 to be con- 
sistent with the high overlap results of [13] . With some 
minor variation of their values, we aim to recover the 
highest overlap possible. Once a good choice of a and 
j3 is found for the 2PN 1.4-1.4 M Q template, we use the 
same values for other templates, but vary the value e (and 
consequently the number of IIR filters in each bank) to 
see the effect on overlap. 



B. Detector Data Simulation 

To test the detection efficiency of the SPUR method 
compared to the frequency domain matched filter, we will 
filter two mock signals, one for which the input data is 
just LIGO-like noise, and the other with the same noise 
plus an inspiral waveform injection scaled to represent a 
source at a chosen effective distance D e ff- 

For this test, we need to construct a finite segment of 
detector data to filter. Because of the IIR filters should in 
principle should be run for an infinite length of the input 
data, we need to run the IIR bank for a finite "warm- 
up" period before the output is consistent with that of 
an IIR filter that has been running for an infinite amount 
of time. In practise, we choose to run each filter for 2 e- 
foldings of time before we accept the output as being 
identical to one which has run for an infinite amount of 
time. Additionally, since each IIR filter in the bank runs 
on a delay, the summed output of all the IIR filters will 
not be produced until after the longest delay time (d max ) 
has passed. The filter that has the longest delay (d max ) 
is also the one that has the longest decay rate 7 max . In 
total, the input data must at least d max + 27"^ in length 
before any output is produced. Hence the length of the 
input data is, 



input 



= d„ 



27n 



N, 



analysis 



(39) 



where -/V ana i ys i s is the length of analysis period, which 
we choose to be 4 seconds. Hence the 4 s SPUR out- 
put will tell us whether there is an injection that ended 



1. Noise generation 

The LIGO-like noise data is produced by creating a 
normally distributed white noise time series of length 
A'input, then colouring it by the theoretical advanced 
LIGO noise spectrum 5 n (/ )[A| We then over- whiten this 
time series using equation (21 ) to produce the waveform- 
free noise input data x: 



'(t). 



(40) 



2. Waveform injection 

We create our waveform injections by first produc- 
ing an inspiral waveform band-limited between 10 and 
1500 Hz. The injection is padded with zeros so that it 
has the length -/Vi nput . The end of the waveform is cho- 
sen so that it finishes somewhere after d m + 27" 1 data 
points. The injection signal is then over whitened using 
equation (21). The over- whitened injection can then be 



placed in the over-whitcned noise signal, 

•Enoise+injection (^) — ^noisc(^) h if )- 



3. Matched filter comparison 



(41) 



As a comparison, we will also perform a frequency do- 
main correlation matched filter. For this process, since 
the input data is already over- whitened, it only needs 
to be cross-correlated with the waveform. Section III Bl 
outlines how this is done. The cosine phase h c (t) gets 
pre-padded with enough zeros to get to length Ai nput . 
This ensures that h c (f) has the same spectral resolution 
as s(f). The matched filter ( 17) produces a time series of 
length. However the first A; nput — Analysis data 



input 



points are erroneous wrap-around caused by the FFT. 
Only the interval [A input - A ana i ysis + 1, A ana i ysis ] is used 
to determine if a waveform is present . 



C. Detection Efficiency 

To test the detection efficiency of the SPUR method 
compared to the traditional matched filter method we 
will construct several receiver operating characteristic 
(ROC) curves for 2PN 1.4-1.4 Mq waveforms injected 
for different effective distances D c g . To create each ROC 
curve, we first find the false alarm rate. The false alarm 
rate is found by realising an Ni nput length LIGO-like noise 
time series, filtering this input data, and analysing the 
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output of the 4 s analysis period (the SNR). We will count 
this realisation as a false positive if at any point within 
the 4 seconds the SNR goes over a given SNR thresh- 
old. Several thresholds will be chosen, giving the false 
positive as a function of threshold. After > 10 6 noise re- 
alisations, the false alarm rate is simply the ratio of total 
number of false positives to number of noise realisations. 
Likewise, to see if the IIR filter doesn't miss too many 
true positives, we inject a 2PN 1.4-1.4 M Q waveform us- 
ing the prescribed method in |III B 2| for a given D e g into 
LIGO-like noise. After filtering, if at any point within 
the analysis period the SNR is above a given threshold, 
this realisation is counted as a true positive. Again, after 
> 10 6 noise realisations, we calculate the detection rate 
as a ratio of the total number of true positives to num- 
ber of realisations. The plot of false alarm rate versus 
detection rate gives the ROC curve. 
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FIG. 4: The overlap between the exact inspiral waveform and 
the approximate inspiral waveform as a function of number 
of damped sinusoids. In general the greater the number of 
sinusoids per waveform, the greater the overlap. However the 
choice of 7;S greatly affects the overlap. 



IV. RESULTS 

A. Inspiral Waveform Overlap 

Starting with the canonical 1.4-1.4 Mq second order 
post-Newtonian binary waveform band limited to be be- 
tween 10 and 1500 Hz we found, using the parameters 
e = 0.04, a = 0.99, f3 — 0.25 in the procedure outlined 
in Section II F that can recover an overlap of 99% using 



687 IIR filters. 

We find that increasing the value of e will in general 
increase the overlap, as the frequency space is more finely 
sampled. However there seems to be a limit, as the damp- 
ing factor 7 causes the adjacent IIR filters to run into 
each other. 

With this choice of a and /3 we are able to recover a 
high overlap for different mass pairs as well. Figure [4] 
shows the overlap as a function of number of IIR filters 
for six different mass pairs. 



- IIR filter output 

- Matched filter output 




FIG. 5: The SNR output of both the SPUR method and a 
traditional matched filter method. The plot is centred on 
t — t c where r c is the time at which the injection ends. From 
the two curves, it is clear that the SPUR method can return 
a very similar SNR to that from the optimal filter. 



B. Ability to Recover SNR 

Figure [5] shows the SNR produced from both the 
matched filter technique and the SPUR method. The 
input time series is constructed following Section |III B| 
The injection of a 2PN 1.4-1.4 Mq waveform scaled for 
an effective distance of 250 Mpc is added to LIGO-like 
noise. The x-axis of the plot is centred about end of the 
injection (t — r c ), which is directly in the middle of the 
analysis period. Around this time, the SNR peaks to 8.2, 
which is near the expected value of 7.9 for an injection at 
this distance. This plot shows that the SPUR method is 
capable of recovering a very similar SNR to the matched 
filter at all times. 



C. Detection Efficiency 

We analysed over 10 6 independent noise realisations, 
for which the waveform had been injected at D c s of 250, 
300, 350, 400 Mpc. We performed both IIR filtering and 
traditional matched filtering. Figure [6] shows that the 
SPUR method recovers most of the same events as the 
traditional matched filter method. At false alarm rates 
of greater than 10 -5 , the SPUR method recovers greater 
than 99% of the injections recovered by the matched fil- 
ter when searching for injections at an effective distance 
of 250 Mpc (SNR~7.9). Even in the worst case, at a false 
alarm rate of 10 -6 , the SPUR method catches 4.5% of in- 
jections scaled at an extreme 400 Mpc (SNR~5), whereas 
the matched filter catches 5% of injections at this scale. 
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FIG. 6: The receiver operating characteristics (ROC) of 
both the IIR filter method and the traditional matched fil- 
ter method. The x-axis shows the false alarm rate, and the 
y-axis the detection rate. A one-to one relationship, which 
is the worst case scenario, is shown by the boundary of the 
shaded area. We show four different ROC curves, where the 
each curve represents the detection rate as a function of false 
alarm rate for waveform injected at effective distances of 250, 
300, 350 and 400 Mpc (SNR ~ 8, 6.6, 5.7 and 5 respectively). 

V. SUMMARY AND DISCUSSION 

The use of a bank of simple IIR filters for each tem- 
plate as opposed to the matched filter method enables 
us get two extra processes for a minimal additional cost. 
The first is that the individual IIR filter outputs can be 
arranged into groups, such that their total summed out- 
put is roughly independent and orthogonal to each other. 
This enables, with minimal extra overhead, the calcu- 
lation of a x 2 distributed statistic, giving a secondary 
method of verification. We will demonstrate this in an 
upcoming paper. The second natural advantage of using 
a parallel bank of single-pole IIR filters is that they can 
easily be executed in parallel using multi-threaded pro- 
cessors, such as graphics processing units (GPUs). In- 
deed, a side study has shown that this is possible [24] . 
This leads to the future possibility that a single personal 
computer may be able to process the detection of GWs. 

A further way to reduce the computation of the IIR 
calculation is to split the incoming data into differently 
down-sampled channels. The output of each IIR filter in 
the bank is the correlation of a fixed frequency sinusoid 
and the incoming data. For the sinusoids that have fre- 
quencies < 124 Hz, the incoming data need only be sam- 
pled at 256 Hz. The current pipeline of LLOID uses a 
similar multi-channel down-sampling in their detection 



pipeline. Their pipeline consists of the integration of 
the open-source real-time multimedia handling software 
gstreamer and the LIGO Algorithm Library (LAL) |14) . 
This software library is an ideal platform to integrate 
the SPUR method. The total computation can also be 
further reduced by sharing IIR filters (via interpolation) 
between different templates [T2] . 

Although the design of the IIR filter so far only ap- 
plies to chirping, post-Newtonian approximation inspi- 
rals, we have performed preliminary tests using more 
complicated combinations of single-pole IIR filters to 
replicate the waveform of an inspiral with spin. If the 
amplitude/frequency beating of a spinning inspiral wave- 
form can be simulated by the linear addition of two dif- 
ferent non-spinning inspirals with different masses, then 
it can be approximated by a linear addition of damped si- 
nusoids. In this case, the SPUR method can produce the 
SNR for the beating waveform. There is also the possi- 
bility of using higher order IIR filters, although designing 
the coefficients can be very difficult. 



VI. CONCLUSION 

We have demonstrated that the through the use of a 
parallel bank of single pole IIR filters, it is possible to 
approximate the SNR derived from the matched filter 
with greater than 99% overlap. The main advantage of 
our SPUR method is that it operates completely in the 
time domain, and in principle it has zero latency (not 
taking into account whitening or computational time). 
The SPUR method recovers most of the injections the 
optimal matched filter recovers. 

We foresee that the use of IIR filters for time do- 
main filtering of Advanced LIGO will be ideal, as the 
waveforms will be much longer. The frequency domain 
matched filter will take more time to calculate GW trig- 
gers, essentially ruling out the possibility of triggering the 
detection the prompt optical emission related to neutron 
star mergers (GRBs). We have shown that the use of 
a parallel bank of IIR filters requires less computational 
cost, with minimal detection rate loss, and most impor- 
tantly can be calculated in the time-domain with near 
zero latency. 
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Appendix A: Noise Spectral Density 

We use an algebraic expression for the noise spectral 
density of Advanced LIGO detectors defined by, 



s h (f) 




(Al) 



where, f = 215Hz and S = 10 
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